Organic magnetoresistance near saturation: mesoscopic effects in small devices 
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In organic light emitting diodes with small area the current may be dominated by a finite num- 
ber, of sites in which the electron-hole recombination occurs. As a result, averaging over the 
hyperfine magnetic fields, bh, that are generated in these sites by the environment nuclei is incom- 
plete. This creates a random (mesoscopic) current component, 51 (B), at field B having relative 
magnitude ^ N~^^^. To quantify the statistical properties of SI{B) we calculate the correlator 
K{B,&B) = {5I{B - ^)SI{B + ^)) for parallel, aB||B, and perpendicular, aB _L B orienta- 
tions of aB. We demonstrate that mesoscopic fluctuations develop at flelds |B| 3> \bh\, where the 
average magnetoresistance is near saturation. These fluctuations originate from the slow beating 
between 5 and To states of the recombining e-h spin pair-partners. We identify the most relevant 
processes responsible for the current fluctuations as due to anomalously slow beatings that develop 
in sparse e-h polaron pairs at sites for which the b,, projections on the external field direction almost 
coincide. 

PACS numbers; 73.50.-h, 75.47.-m 
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Introduction. In the field of 'Dynamic Spin Chemistry', 
a mechanism by which the recombination rate of radical 
pairs is sensitive to a weak magnetic field, B, was estab- 
lished more than four decades ago; see, e.g., the reviews 
m Ref. m This mechanism relies on the hyperfine inter- 
action of the spin-1/2 pair partners with their respective 
nuclear spin environments, where the hyperfine field, 
generated by the nuclei is responsible for the radical spins 
dynamics in zero field. In this process if at time t = the 
radical pair spin state is, e.g., in a singlet configuration, 
5, then at finite t it will acquire a triplet (T) component 
with probability, Psrit). If recombination is allowed only 
from S, then Psr(0 dynamic evolution affects the recom- 
bination rate. Clearly, Psrit) depends on B and this sets 
a small scale, |B| ^ jb^l that may influence the radical 
pair recombination rate. 

An important advance in the quantitative description 
of PsT (i) was made by Schulten and Wolynes^. They no- 
ticed that, due to the large number of nuclei surrounding 




FIG. 1: The dependence I{B) of the device current on the ap- 
plied magnetic field is shown schematically in the strong-field 
limit B ^ bo. Enlargement illustrates mesoscopic fiuctuations 
emerging in a small sample. Two insets are the correlators of 
the mesoscopic fluctuations for aB |[ B and aB _L B plotted 
from Eqs. (|f 4|l and (|16p . respectively . 
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FIG. 2: In strongly inhomogeneous device with W ^ L the 
current passage is dominated by the most conductive chan- 
nels, / = X^n-'^"- Each current component is limited by the 
most resistive junction, illustrated schematically. The cur- 
rent through this junction is sensitive to the spin dynamics of 
the constituting PP. "Slow" pairs, shown in the enlargement, 
are those in which the ^-projections of their hyperfine fields 
coincide accidentally. 



each radical pair, and slow dynamics of the hyperfine 
field, the random distribution may be modeled by a 
Gaussian. Under these conditions the multiplicity of the 
nuclear spin configurations may be characterized by a 
single number - namely the width of this distribution, 
bo- 

The dependence of Psrit) on B is at the core of or- 
ganic magnetoresistance (OMAR), which has recently at- 
tracted a lot of attention--!^. This is because the cur- 
rent, / in a biased organic diode involves recombination 
of the injected e-h polarons (forming polaron pairs, PP); 
whereas the processes of populating and depopulating of 
traps are not sensitive to spin dynamics. 

The theory of OMAR is conceptually harder than that 
of spin-magneto-chemistry-'^ for two reasons. Firstly, in 
OMAR the complex dynamics of all four PP spin states, 
5, To, T+, and T_ needs be incorporated into the cal- 
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culation of the dc current that is influenced by the PP. 
Secondly, each PP is sensitive to the other PPs if they 
belong to the same current path. Finally, averaging over 
the nuclear environment should be carried out only at 
the last step. To bypass these complications several sim- 
plifying assumptions concerning both the spin dynamics 
and current passage scenario were adopted in previous 
theoretical calculations of the I{B) response^iii"— . 

In contrast, in the present paper we do not focus on 
the entire I{B) response, but rather on the strong B do- 
main, where the OMAR response is close to saturation, 
see Fig. [TJ Our motivation is twofold. Firstly, theory 
allows a dramatic simplification in this i3-domain, since 
the spin dynamics that is relevant to OMAR involves 
only the PP S and Tq states. However even in this B- 
domain the OMAR underlying physics is not trivial if 
the hyperfine field is sufficiently strong; namely when 
boT ^ 1, which corresponds to the regime of "slow" 
hoppingi^iii. Here r is a characteristic recombination 
time. Experimentally^^, in organic semiconductors 60 is 

1 mT, whereas r ~ 1 — 10 ^s, so that this parameter is 
^ 10"^. We show that at large hgT the spin dynamics is 
not "frozen" as B exceeds 60, but persists in a paramet- 
rically broad interval, b^r » i? » 6o- Our second and 
central motivation for considering strong fields is that we 
predict the occurrence of mesoscopic properties in this B- 
domain that would form in small devices that are based 
on strongly disordered organic active layers. Specifically 
we predict reproducible random fiuctuations in the I{B) 
response upon sweeping B (see Fig. [T|) , which reflect the 
"individuality" of the nuclear environments associated 
with the relevant recombination centers in the organioiS,. 

More quantitatively, if the number, N of current paths 
that contribute to I{B) is finite, then the statistical av- 
eraging over 6h is incomplete. The relative fluctuation 
^iiyl iV^^/^ while small, can be still experimentally 
obtained because of the high accuracy with which current 
can be measured. In the fleld of 'dynamic spin chemistry', 
mesoscopic fluctuations cannot occur since the number of 
radical pairs that contribute to the observable character- 
istics is huge. 

Obviously, the necessary condition to observe meso- 
scopic fluctuations in the I{B) response of organic de- 
vices is slow nuclear spin dynamics, which should allow 
to obtain I{B) before the nuclear configuration changes. 
This is realistic, since the characteristic time of nuclear 
spin change is ~ 10'* times longer than the PP recombi- 
nation time, r due to the weakness of the nuclear dipole- 
dipole interaction in organic semiconductors. 

We start with a detailed account of the PP spin dy- 
namics and recombination in the strong i3-domain, which 
we then use to calculate I{B) near saturation. 
S — To beating in the strong-field limit. For an isolated 
PP the spin Hamiltonian H = fli- Si + n2- S2 describes 
the precession of the PP spins 5'i, ^'2 in the fields Oi = 
B+bi and = B+b2, respectively. If at i = the PP is 
in the singlet state, then the probability, Pss{t) to find it 
in the singlet state at finite t oscillates with time. Pss (t) 



oscillations contain two frequencies: A = — |r22| and 
S — |r2i|-)-|r22|. The advantage in considering the strong- 
field limit is that since \ w \^2\, the frequencies A 
and S are very different from each other, so that the 
spin dynamics decouples into distinct 'slow' and 'fast' 
modes. Moreover, the slow mode involves predominantly 
S and Tq states, while the admixture of T_^_ and T_ states 
to this mode is relatively weak (of the order of b^/B'^). 
The fast mode E has frequency w 2B and describes the 
oscillations between 5* and T^, T_. But the admixture 
of S to this mode is also suppressed as b^/B"^ in the 
strong-field limit. We thus conclude that, with accuracy 
bfjB^, Pss{t) dynamics simplifies in the strong-field limit 
to Pss{t) = cos^ At; namely the 'beating' between S and 
To states. Similarly, if in the strong-field limit the PP is 
initially in the Tq state then the probability to find it in 
the S state at time t is sin^ At. 

Recombination in the presence of S — Tg beating. We now 
assume that the PP is still isolated from the 'leads', but 
can recombine from S to the ground state, G. A crucial 
question for OMAR is: what are the waiting times {t)s, 
{t)TQ for the recombination, if the system is initially in 
S and To, respectively. The simplified spin dynamics in 
the strong-field limit allows us to address this question 
analytically. 

Upon restricting the basis to S, Tg and the ground 
state, we have 9 relevant elements of the density matrix 
for solving the Liouville-Lindblad equations of motion: 
p = —i[H,p] + L{p), where the operator L{p) describes 
the recombination. To find, e.g. {t)s the system should 
be solved with the initial conditions p{0) — \S) {S\. Sub- 
sequently {t)s is found from the formula 

{t)s^ dtt-^=J^ dt{pss{t)+PT,TM- 

Similarly (i)r(, is obtained from Eq. ([1]) upon solving 
the equations of motion with initial conditions p(0) = 
I To) (To|. These calculations yield 



{t)s^r, {t)^^=r+^. (2) 

For a typical PP we have A = |r2i| — |J72| bg. Eq. 
([U suggests that {t)s ~ (i)To ~ This is a natural 
result since recombination is preceded by many beatings 
between S and Tq states; therefore the recombination 
time does not depend on the initial PP state. The most 
striking consequence of Eq. ^ is that for sparse PP 
for which A is accidentally smaller than we have 
{t)s ^ {t)Tg ~ 2A^r • '^^^^ suggests that the smaller is 
A, the longer the pair stays "trapped" in Tg. Note that in 
the course of beating without possibility of recombination., 
such PP would cross from Tg to S after much shorter time 
A~^ <C {t)To- We can trace the origin of the "trapping" 
described by Eq. ([2]) to the complex eigenmodes of the 
system that consists of singlet and triplet components 
being mixed by the hyperfine field. This system may be 
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FIG. 3: Schematic representation of the PP population dy- 
namics in the strong-field domain. For variants of cycle I (IV) 
the pair is assembled and, subsequently, disassembled in T+ 
(T_) state. For variants II (III) the pair is assembled in S 
(To) state in which it undergoes slow dynamics prior to dis- 
assembly. In the course of the slow dynamics the pair can 
recombine; recombination is possible only from S. Since the 
transport is unidirectional, current is passed through a junc- 
tion upon completion of each cycle variant. 



described by the 2x2 non-hermitian matrix: 
1^) \T) 
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(3) 



where the nondiagonal elements describe the mixing, 
while — z/r describes recombination from S to G. The 
eigenvalues of this matrix are 
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In the limit A <C r ^ we have Ai 



while A, 



jrA^. We see that Aa is anomalously small, and the result 
Eq. for A <C can be interpreted as {i)To ^ x;- 
We note in passing that the emergence of slow mode, 
A2 in a compound system with anomalously close levels 
was previously found in Refs. [l9ll20l in connection with 
resonant tunneling through pairs of localized states. 

In the remainder of the paper we demonstrate that it 
is the sparse pairs with A < r""'^ that are responsible for 
/(B) response in the strong-field limit, as well as for the 
mesoscopic fluctuations in /(-B) response. 
Transport model. We assume that the organic active layer 
in the device is strongly inhomogeneous and its width, 
W is much larger than the distance, L between the elec- 
trodes, see Fig. [2j Under these conditions there are 
^ 1 parallel conducting channels. Moreover due to 
the film inhomogeneity, the current through each chan- 
nel is limited by a single, most resistive junction. The 
stronger is the inhomogeneity, the more realistic is the 
proposed model. The net current through the sample is 
the sum, / = J2n=i ^n, of the currents in each channel. 
Each junction, n, can be viewed as a pair of sites coupled 
to the nuclei environment. In the course of the current /„ 



through the junction, the pair of sites first gets occupied, 
is then emptied, and so on. In other words, the current 
passage can be viewed as a sequence of cycles, see Fig. 
[21 Each cycle consists of two steps, namely assembly of a 
pair on neighboring sites and disappearance of the pair 
through either dissociation or recombination, see Fig. [3l 
At this point we emphasize that it is the recombination 
stages of the cycles II and III (defined in Fig. [3]) that are 
described by Eq. ([2]), and are thus sensitive to B. 

From Fig. |3]we may write the average duration (i„) of 
the cycle as follows 

(t.) = 2x i(2r„) + i((t),+T^) + i((i),„+r„). (5) 

The first term in Eq. ([S|) originates from the variants I 
and IV of the current cycle when the pair is assembled, 
and subsequently disassembled in and T_ states, re- 
spectively. Then the cycle lasts for time 2to- The last 
two terms in Eq. ([5]) describe the current cycle variants 
II and III, in which the pair is assembled in either S' or Tq. 
Eq. ([5]) takes into account that realization of each current 
cycle has equal probability of ^ . For simplicity we assume 
that processes involving leads, namely, assembly and dis- 
sociation, take equal time, Td- The current through the 
junction, which is the inverse cycle duration, can be then 
cast in the form 
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8C^T-dIJB), SI„{B) 



8CV 



A„{Bf + e' 
(6) 

where we used Eq. ([2]) and introduced the characteristic 
frequency 



c 



2[t{3to+t)] 



1/2- 



(7) 



We emphasize that the correction, SI„{B), in Eq. ^ 
originates from S — beating. For a typical nuclear 
environment we have A„ ^ bg, so that the relative mag- 

nitude of this correction is ~ W <C 1. However, on av- 

° . 

erage, this term is much bigger, since it is dominated by 
sparse configurations with anomalously small A„ ^ C,. 
This is because, while the portion of these sparse con- 
figurations is small, ^ C/^o, the (5/„ value for these 
configurations exceeds the typical 51^ by a large factor 
^. The latter observation allows us to find the aver- 
age {61(B)), since in calculating this average, the factor 
(A2-HC^)-i can be replaced by fS{\B + bi\-\B + b2\). 
Averaging of the (5-function over the components bi of 
the random hyperfine fields, each of which is distributed 
as 'P{bi) = -y=5^ ^^p{~bf /bg), can be performed analyti- 
cally. In the strong-field limit this averaging yields 



I{^) {IjB)) 

/((X)) 



27r boC 
8 B2 



^ f- (8) 



Note that for a typical PP the argument of the (S-function 
in the strong-field limit contains a correction ^ 1/ B. The 
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overall 1/_B^ dependence of Eq. ([8]), illustrated in Fig. 
[T] describes the fact that these corrections have random 
signs, and cancel out on average. The factor Q/ha in Eq. 
([5]) reflects the portion of the sparse hyperfine configu- 
rations with slow spin dynamics. This slow dynamics is 
so explicitly reflected in OMAR via Eq. ([8|) due to the 
extreme simplicity of the model adopted. In a more so- 
phisticated treatmentii, which incorporates spin dynam- 
ics in calculating OMAR, the dominance of pairs with 
slow dynamics is not so explicit as in Eq. ([8]), but is indi- 
rectly present in the slow-hopping limit. Eq. ([8|) suggests 
that the spin dynamics enters into the net current only 
crudely, namely, via the ratio, C/^o) of the "slow" pairs 
involved. This dynamics manifests itself much more del- 
icately in the mesoscopic properties of OMAR, as seen 
below. 

Mesoscopic fluctuations. If a given PP pair contributes 
to the correction 6I„{B) in Eq. (jG)), then the S — Tq split- 
ting, A„(B) for this pair is ^ C- This suggests that, upon 
changing B by a small &B, the condition A„(B) ^ C for 
this pair is violated, while it becomes satisfied for dif- 
ferent pairs. Such "switching" of pairs contributing to 
the correction, SI{B) gives rise to the mesoscopic fiuc- 
tuations of the current, which we may quantify by the 
correlator 

KiB,.B) ^ (^SliB- ^)6I{B + 

(9) 

We consider two cases. In the first case, ais || B, the two 
magnetic fields are coUinear. In the second case, ai? _L B, 
the two magnetic fields have the same magnitude but are 
rotated through an angle ip with respect to each other. 

For calculating the correlator in Eq. ([9|) we recast the 
factor (A^ + C^)^^ in <5/„ as a Fourier transform 
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A2 + ^2 -^y_^^^exp(-|,s|+*— ). (10) 
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By virtue of this transformation, the beating frequency 
A, which depends on the hyperfine fields, appears in the 
exponent of the integrand. Next we take advantage of 
the fact that the beating frequency, A in the strong-field 
limit can be expanded as 



A^{b,~b^)-n+ — , (11) 



where n is the unit vector in the direction of B. Now, 
since A contains only linear and quadratic terms in bj, 
the averaging of the exponential factor can be performed 
explicitly using the properties 



(12) 



For the parallel case, this averaging yields 



{A{B + ^)2 + e){A{B - ^)2 + C 



4^4 



dsi / ds 



exp (^-(|si| + Is^l) - ^ (si - Sa)^) 




(13) 



1 



/ Si 



The integration over the difference s^ — Sa can be per- 
formed explicitly using the fact that ( <^ bg ^ B. Upon 
this integration, the average Eq. (1131) can be presented in 
$11 (^), where the dimensionless function 



the form 
$11 is defined as 



dx 



1 -|- Z^X'^ 



(14) 



and = AB^(/bg. The argument of the function $ im- 
poses a characteristic period of mesoscopic fluctuations: 
^ ~ Our consideration applies for SB <^ B. This 

condition suggests that for measuring the fluctuations 
one must work in the domain bg ^ B ^ &o/C- The cor- 
relator Eq. ([H]) is plotted in Fig. [T] For smah SB < B^ 
it behaves as 1 — (^)^, and falls off slowly, as for 
SB > B,. 

For the perpendicular case, we can simplify A as 
A « {b^ — b^) ■ n. This is because the B-dependence 
of A enters via the orientation, n. Performing the same 
decoupling (Eq. (fTO| ) as for the parallel case, instead of 
the double integral in Eq. (IT^ we get now 



dsi 



ds2 



- 1 '2 1 exp <{- ^ (s^ s2 _ cos vs) 



(15) 

We again see that by virtue of the relation ho ^ C, the 
difference (si — Sg) ~ ^ is small. This allows us to inte- 
grate over Si — Sa and reduce Eq. ([15]) to ^^|^ $_l(^), 
where ip^ — ^r-^ ^ 1 and the function '^a_{z) is defined 



as 



poo 

/ dx exp[- 

"'0 



2 21 

■X — Z X 



(16) 



The correlator is plotted in Fig. [TJ At 93 ^ (p^ , it falls off 
as ^Pa/'P>- In general, the correlator, Eq. (fTS)) is a periodic 
function of (/?; had we not used the small-(^ expansion it 
would go through a minimum at ip = Tr/2 and "revive" 
at if = TT. 

Our results related to mesoscopics can be summarized 
in the following expression 

(SI{B-^)SI{B + ^)) ^ f*,, (^) ,^B\\B 

^(oo)2 mbo 1 $, (^] .B±B 



(17) 

Concluding remarks. By choosing a simple transport 
model for an organic semiconductor device, and making 
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the strong assumption that recombination proceeds ex- 
clusively from the singlet state, we were able to capture 
analytically the OMAR response in the domain B ^ 60, 
and predict the characteristic magnitude and period of 
mesoscopic fluctuations. Our theory is based on an im- 
portant finding that in this _B-domain there exists a 
strong separation between slow and fast components of 
the PP spin-dynamics. As a result of this separation, the 
S-Tq beating becomes decoupled, which, in turn, leads 
to a dramatic recombination slow down which originates 
from PP "trapping" in the Tq state. Since this underly- 
ing physics is so general, any transport model in a small 
device with few junctions should exhibit mesoscopic fluc- 
tuations. What is really required for mesoscopic features 
to emerge in I{B) is that the transport is in the regime of 
"slow-hopping" , namely bo 3> C; when the local current 
is dominated by sparse configurations of nuclear environ- 
ments. In the opposite limit, namely 60 ^ C: can 
see from Eqs. © and ([TT|) that OMAR response in the 

strong- field domain has the form ^ How- 

ever Eqs. (fT3|) and (fT5| do not yield mesoscopic fluctu- 
ations in this case. This could be expected on general 
grounds, since bg <ti C means that separation of the Zee- 



man levels is much smaller than their broadening. 

The simplified transport model described in Fig. [2]en- 
codes the same picture of transport put forward in Ref. 
I12I . It contains recombination and bypassing, if recombi- 
nation takes too long. However, the "probabilistic" de- 
scription adopted in Ref. [H excludes the current sensi- 
tivity to the spin dynamics which, as we showed here, is a 
necessary precondition for obtaining mesoscopic physics. 

Regarding experimental verification of the predicted 
mesoscopic fluctuations, we note that ITO electrodes ex- 
hibit sharp pillars^^ with areal density of ~ l/im~^. 
These pillars may cause high-conductivity channels. 
Then a small OLED with area of ^ 10~^ cm^ will show 
mesoscopic fluctuations of y" ^ 10^^. 

As a final remark, note that the slow beating is char- 
acterized by the time scale r, the time required for re- 
combination. However, transport characteristics contain 
1/C given by Eq. [3 which contains both t and To- The 
latter being the measure of the PP coupling to the envi- 
ronment (electrodes). For ^ t, ( is the geometrical 
mean, C ~ (tto)^^/^. 
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